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Abstract: We consider the resummation of soft gluon emission for squark-antisquark 
pair-production at the LHC at next-to-next-to-leading-logarithmic (NNLL) accuracy in 
the framework of the minimal super symmetric standard model. We present the analytical 
ingredients needed for the calculation and provide numerical predictions for the LHC at 
centre-of-mass energies of 7 and 14 TeV. We find a significant reduction in the scale uncer- 
tainty and a considerable increase in the prediction of the total cross section. Compared 
to the next-to-leading order prediction, the corrections increase the cross section by up to 
30% for 1.5 TeV squarks at a centre-of-mass energy of 7 TeV. 
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1. Introduction 

The search for supersymmetry (SUSY) [1,2] is a central part of the experimental program at 
the Large Hadron Collider (LHC). Models of weak-scale SUSY provide a promising solution 
to the hierarchy problem of the Standard Model (SM) and predict new supersymmetric 
particles (sparticles) with masses in the TeV range. If they exist, the coloured sparticles, 
squarks (q) and gluinos (<?), would be produced copiously in hadronic collisions and thus 
offer the strongest sensitivity for supersymmetry searches at the LHC. In the context of the 
minimal supersymmetric extension of the Standard Model (MSSM) [3,4], these particles 
are produced in pairs due to imparity conservation. 

Searches for squarks and gluinos at the proton-proton collider LHC, which has been 
operating at y/S = 7 TeV in 2010 and 2011, have placed lower limits on squark and gluino 
masses around 1 TeV [5,6]. Once the LHC reaches its design energy of y/S = 14 TeV, 
SUSY particles with masses up to 3 TeV can be probed [7,8]. 

Accurate theoretical predictions for inclusive squark and gluino cross sections are 
needed both to set exclusion limits and, in case SUSY is discovered, to determine SUSY 
particle masses and properties [9, 10]. The inclusion of higher-order SUSY-QCD correc- 
tions significantly reduces the renormalization- and factorization-scale dependence of the 
predictions. In general, the corrections also increase the size of the cross section with re- 
spect to the leading-order prediction [11,12] if the renormalization and factorization scales 
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are chosen close to the average mass of the produced SUSY particles. Consequently, the 
SUSY-QCD corrections have a substantial impact on the determination of mass exclusion 
limits and would lead to a significant reduction of the uncertainties on SUSY mass or pa- 
rameter values in the case of discovery. The squark-antisquark production processes have 
been known for quite some time at next-to-leading order (NLO) in SUSY-QCD [13, 14]. 
Electroweak corrections to the 0{a%) tree- level processes [15] and the electroweak Born 
production channels of 0(aa s ) and 0{a 2 ) [16] are in general significant for the pair pro- 
duction of SU(2)-doublet squarks qi and at large invariant masses, but they are moderate 
for inclusive cross sections. 

A significant part of the NLO QCD corrections can be attributed to the threshold 
region, where the partonic centre-of-mass energy is close to the kinematic production 
threshold. In this region the NLO corrections are dominated by soft-gluon emission off 
the coloured particles in the initial and final state and by the Coulomb corrections due 
to the exchange of gluons between the slowly moving massive sparticles in the final state. 
The soft-gluon corrections can be taken into account to all orders in perturbation theory 
by means of threshold resummation techniques [17,18]. The Coulomb corrections can be 
summed to all orders by either using a Sommerfeld factor [19-22] or by employing the frame- 
work of non-relativistic QCD, where bound-state effects can be included as well [22-27]. 
In addition, a formalism has been developed in the framework of effective field theories 
that allows for the combined resummation of soft and Coulomb gluons in the production 
of coloured sparticles [28,29]. 

Threshold resummation has been performed for all MSSM squark and gluino pro- 
duction processes at next-to- leading-logarithmic (NLL) accuracy [30-34]. For squark- 
antisquark production, in addition to soft-gluon resummation, the Coulomb corrections 
have been resummed both by using a Sommerfeld factor [31] and by employing the frame- 
work of effective field theories [29]. Furthermore, the dominant next-to-next-to-leading or- 
der (NNLO) corrections, including those coming from the resummed cross section at next- 
to-next-to-leading-logarithmic (NNLL) level, have been calculated for squark-antisquark 
pair-production [35]. 

In this paper we consider threshold resummation at NNLL accuracy for squark-anti- 
squark pair production at the LHC. Compared to the NLL calculation the new ingredients 
are the one-loop matching coefficients, which contain the NLO cross section near threshold, 
and the two-loop soft anomalous dimensions. Studies for the pair production of top quarks 
suggest that the effect of the matching coefficients can be significant [36] and that NNLL 
resummation can reduce the scale dependence considerably [35,37-40]. We will discuss the 
impact of the corrections and provide an estimate of the theoretical uncertainty due to 
scale variation. In addition we will study the impact of the Coulomb gluons on the cross 
section. We exclude top squarks from the final state in view of potentially large mixing 
effects and mass splitting in the stop sector [41]. The other squarks are considered as being 
mass degenerate and all flavours and chiralities are summed over. 

The structure of the paper is as follows. In section || we discuss the NNLL resumma- 
tion for squark-antisquark pair-production. In section [| we present the calculation of the 
hard matching coefficients required for the NNLL resummation. The numerical results are 
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presented in section |j. We show predictions for the LHC with centre-of-mass energies of 
y/S = 7 TeV and y/S = 14 TeV. We will conclude in section |. 



2. Threshold resummation at NNLL 

In this section we briefly review the formalism of threshold resummation for the production 
of a squark-antisquark pair. The inclusive hadroproduction cross section a hlh2 ^g= can be 
written in terms of its partonic version a^^— as 

<y hl h 2 ^{pA m2 }) = J dxidx2d P 5 {p~^^) 

x fi/hi( x l, P 2 ) fj/h 2 ( x 2, P 2 ) ^^{^{m 2 }, P 2 ) , (2-1) 



where {m 2 } denotes all masses entering the calculations, i,j are the initial parton flavours, 
fi/hx and fj/h 2 the parton distribution functions and p is the common factorization and 
renormalization scale. The hadronic threshold for inclusive production of two final-state 
squarks with mass rriq corresponds to a hadronic centre-of-mass energy squared that is 
equal to S = 4m?. Thus we define the threshold variable p, measuring the distance from 
threshold in terms of energy fraction, as 

4m? 

The partonic equivalent of this threshold variable is defined as p = p/(x\X2), where xi 5 2 
are the momentum fractions of the partons. 

In the threshold region, the dominant contributions to the higher-order QCD correc- 
tions due to soft-gluon emission have the general form 

4m? 

a? log m fi 2 , m < 2n with f = 1 - p = 1 , (2.2) 

where s = X1X2S is the partonic centre-of-mass energy squared. The resummation of the 
soft-gluon contributions is performed after taking a Mellin transform (indicated by a tilde) 
of the cross section, 

^ hl h 2 ^( N A m2 }) = I d P P N ~^ e^h^fqipA™ 2 }) 
Jo 

= /vh 1 (A r + l,M 2 )/ i /h 2 (A r + l,^ 2 )^^(^W 2 }^ 2 ) • ( 2 - 3 ) 

The logarithmically enhanced terms are then of the form a™ log m N, m < 2n, with the 
threshold limit (3 — > corresponding to N — > 00. The all-order summation of such loga- 
rithmic terms is a consequence of the near-threshold factorization of the cross sections into 
functions that each capture the contributions of classes of radiation effects: hard, collinear 
(including soft-collinear) , and wide-angle soft radiation [17,18,36,42-44]. Near threshold 
the resummed partonic cross section has the form: 
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x Ai(JV + 1, Q 2 ,/x 2 ) A,-(iV + 1, Q 2 , » 2 ) Ag^ v (Q/(iV/x), /i 2 ) , (2.4) 

where we have introduced the hard scale Q 2 = 4m|. Before commenting on the different 
functions in this equation, we recall that soft radiation is coherently sensitive to the colour 
structure of the hard process from which it is emitted [36,42-46]. At threshold, the resulting 
colour matrices become diagonal to all orders by performing the calculation in an s-channel 
colour basis [28,30,31]. The different contributions then correspond to different irreducible 
representations I. For the qq — > qq process, the s-channel basis consists of a singlet 1 
and an octet 8 representation, while for the gg —> qq process it contains a singlet 1, 
an antisymmetric octet 8a and a symmetric octet 8s representation as presented in e.g. 
Ref. [28,30,31]. 

In Eq. fl2.4|), a^'_^.= j are the colour-decomposed leading-order (LO) cross sections in 
Mellin-moment space. The functions Aj and Aj sum the effects of the (soft-)collinear 

(si 

radiation from the incoming partons, while the function A . _= , describes the wide-angle 

& 1 ' i]^qq,I & 

soft radiation. Schematically the exponentiation of soft-gluon radiation takes the form 



[17,18] 

Lgi(a s L) + g 2 {a s L) + a s g 3 (a s L) + . . . . (2.5) 



AiA,A (s) _ B , = exp 

1 J ij^tqq,! ^ 



This exponent captures all dependence on the large logarithm L = logiV. Keeping only 



the <7i term in Eq. (2.5) constitutes the leading logarithmic (LL) approximation. Including 
also the g 2 term is called the NLL approximation. For the NNLL approximation also the 53 
term needs to be taken into account. Explicit expressions for the 53 term and its ingredients 
are given in Refs. [47-49] 1 and are listed in appendix^. 

We also need the matching coefficients C, which contain the Mellin moments of 
the higher-order contributions without the log(iV) terms. To NNLL accuracy, this non- 
logarithmic part of the higher-order cross section near threshold factorizes into a part that 
contains the leading Coulomb correction C Coul '^ and a part that contains the NLO hard 
matching coefficients [29]: 



C NNLL = (l + ^C Co ^(N,{m 2 },n 2 )) (l + ^C«({m 2 },/i 2 )) (2.6) 



The calculation of the NLO hard matching coefficients and the Coulomb contribu- 
tion C Coul 'W for the squark-antisquark production processes will be discussed in detail in 
section ||. 

Having constructed the NNLL cross-section in Mellin-moment space, the inverse Mellin 
transform has to be performed in order to recover the hadronic cross section a hlh2 ^q=. In 
order to retain the information contained in the complete NLO cross sections [14], the 
NLO and NNLL results are combined through a matching procedure that avoids double 



1 One has to correct for an extra minus sign in front of all Dqq terms in Eq. (A9) of [49] 
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counting of the NLO terms: 

(NLO+NNLL matched) / r 2i 2\ (NLO) / r 2t 2\ fr, ~\ 



■ ■ J CT 



We adopt the "minimal prescription" of Ref. [50] for the contour CT of the inverse Mellin 
transform in Eq. ( |2.7| ). In order to use standard parametrizations of parton distribution 
functions in x-space we employ the method introduced in Ref. [51]. 

3. Calculation of the matching coefficients 

In this section we will discuss the calculation of the matching coefficients C at one loop. 
As discussed in equation ( |2.6| ), the NNLL matching coefficient C NNLL factorizes into the 
Coulomb contribution and the hard matching coefficient. For NNLL resummation, both 
terms are needed up to NLO accuracy. 

The terms in the NLO cross section which give rise to the Coulomb corrections C Conl '^ 
in iV-space do not have the usual phase-space suppression oc /3, in view of the Coulombic 
1//3 enhancement factor. After performing an expansion of the NLO cross section in (3, 
the hard matching coefficients are determined by the terms in the NLO cross section 
that are proportional to (3, /31og(/3) and /31og 2 (/3). Terms that contain higher powers of 
/3 are suppressed by powers of 1/N in Mellin-moment space and do not contribute to the 
matching coefficient C. In contrast to the case of top-pair production in Ref. [52], there 
is no full analytic result for the real corrections to squark-antisquark production, so we 
cannot take the explicit threshold limit. For the virtual corrections, which also contain the 
Coulomb contribution, we will use the full analytic expressions, but for the real corrections 
we need a different approach. 

To obtain the virtual corrections for squark-antisquark production we start from the 
full analytic calculation as presented in Ref. [14]. As described in detail in Ref. [14], 
the QCD coupling a s and the parton distribution functions at NLO are defined in the 
MS scheme with five active flavours, with a correction for the SUSY breaking in the MS 
scheme. The masses of squarks and gluinos are renormalized in the on-shell scheme, and 
the top quark and the SUSY particles are decoupled from the running of a s . 

To obtain the virtual part of the hard matching coefficients, we first need to colour- 
decompose the result and then expand it in /3. For the first step we only need the colour 
decomposition of the LO matrix element. Due to the orthogonality of the s-channel colour 
basis, the full matrix element squared is then automatically colour-decomposed: 

|A*Inlo,i = 2Re(A^NLoA^ 0i/ ). 

We are now left with an expression in terms of masses, Mandelstam variables and scalar 
integrals. Since we need the cross section to C(/3), we have to expand |A4| 2 to zeroth order 
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in (3. For the squark-antisquark production processes the factors that multiply the integrals 
do not contain negative powers of (3, so we do not have to expand the scalar integrals beyond 
zeroth order in (5. 

The number of integrals that need to be expanded can be reduced. By using the fact 
that the two outgoing momenta are equal at threshold, we can reduce some of the three- 
and four-point integrals to two- and three-point integrals respectively. This procedure can 
be used only for integrals that contain both outgoing momenta. The result of the remaining 
integrals is explicitly expanded to zeroth order in /3. 

Special attention has to be paid to the Coulomb integrals. First, in order to calcu- 
late the Coulomb corrections C Conl '^ in iV-space, we need to know the Coulomb part of 
the NLO correction in /3-space, corresponding to the leading terms in f3 of the Coulomb 
integrals. These leading terms are given by [19-22]: 



Coul,(l) ^ _"s Tf_ (o) 

with k colour coefficients that depend on the process and the dimension of the representa- 
tion. For the (/(/-initiated process they are given by [31]: 

while for the (/(/-initiated process they are: 



3' -gg^qq^A g anc ^ K gg-+qq,8 s g 



K gg^qq,l ~ o> K gg^qq,S A ~ c an( ^ ^gg-^qq^s 



The Mellin transform a Coul '^ of Eq. (3.1) is presented in appendix^. The function C Conl '^ 
can be obtained by dividing a Conl '^ by the Mellin transform of the LO cross section, which 
can be found in Ref. [31]. 

Secondly, the next term in the ^-expansion of the Coulomb integrals contributes to 
the hard matching coefficients. Due to their 1/(3 behaviour Coulomb integrals cannot be 
reduced to lower-point integrals, so they need to be expanded explicitly. 

To obtain the integrated real corrections at threshold, the key observation is that 
they are formally phase-space suppressed near threshold unless the integrand compensates 
this suppression. Therefore we can construct the real corrections at threshold from the 
singular behaviour of the matrix element squared, which can be obtained using dipole 
subtraction [53,54]. We will briefly review the procedure of dipole subtraction and specify 
how only the singular contributions survive in the threshold limit. 

Dipole subtraction makes use of the fact that the cross section can be split into three 
parts: a part with three-particle kinematics , one with two-particle kinematics , 
and a collinear counterterm a c that is needed for removing the initial-state collinear sin- 
gularities. These parts are well-defined in n = 4 — 2e dimensions, but their constituents 
diverge for e — > 0. With the aid of an auxiliary cross section a A , which captures all singu- 
lar behaviour, all parts are made finite and integrable in four space-time dimensions. This 
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auxiliary cross section is subtracted from the real corrections cx R at the integrand level to 
obtain and added to the virtual corrections <r v , which defines 

= J [da R - da A ] e=Q + J [da v + J da A ] e=Q + a c = a™ + a™ + a c 

We will first argue that we can neglect . Compared to the case of two-parton kinematics, 
the phase space of is limited by the energy of the third, radiated massless particle. 
Near the two-particle threshold, the maximum energy of the radiated particle, and thus 
the available phase space, equals E max = ^/s — 2rriq oc /3 2 . Since after subtracting o A no 
divergences are left in the integrand of a* 3 }, the leading contribution of <x< 3 > is at most 
proportional to f3 2 and can thus be neglected. This leaves us with: 

^.NLO.thr _ a {2},th^ _|_ ^C.thr _ ^V.thr _|_ ^.C.thr _|_ ^.A.thr 

so at threshold the real radiation is indeed completely specified by the singular behaviour 
contained in a A . In Ref. [54] the general form of a A is determined by summing over dipoles 
that correspond to pairs of ordered partons. These dipoles describe the soft and collinear 
radiation and reproduce the matrix element squared in the soft and collinear limits. To 
obtain the cross section, the dipole functions need to be integrated over phase space and in 
particular over the momentum fraction x that is left after radiation. In the threshold limit 
the available phase space sets the lower bound of the ^-integral to 1 — /3 2 , while the upper 
bound equals 1. Therefore we cannot get a result of 0{j3) unless the integrand diverges at 
x = 1, which is the case only for soft-gluon radiation. As a result we only need to take into 
account the dipoles that describe gluon radiation. 

Special attention has to be paid to the massive final-state dipole function. In Ref. [54] 
this dipole function has been rewritten in order to simplify the integration. Unfortunately 
this results in a deformation of the phase space integration which changes exactly the finite 
terms that we are looking for. Therefore the expression given in Eq. (5.16) of Ref. [54] 
cannot be used for our calculation and we have to use the original dipole function instead. 
A more detailed argument can be found in appendix |C]. The divergent part of the dipole 
function is completely determined by the soft limit. Using the eikonal approximation we 
obtain the final-state dipole function that correctly reproduces the soft limit at threshold. 
Its behaviour is given by: 

where p g is the gluon momentum and the sums run over final-state particles with momenta 
Pjj, masses m^y and colour charge operators that are defined in Ref. [54]. This ex- 
pression vanishes at threshold, so the final-state dipoles do not contribute in the threshold 
limit. For the other dipole functions and the collinear counterterm we can use the equations 
in [54] and take the threshold limit. 

After combining the real and the virtual corrections the hard matching coefficients can 
be obtained by taking the Mellin transform and omitting the Coulomb corrections and the 
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log(iV) terms. The complete expressions for the hard matching coefficients of the squark- 
antisquark production processes can be found in appendix Their behaviour for varying 
gluino mass is shown in Fig. |[ For the gg — > qq process the antisymmetric octet 8a 




3000 



(a) m i (b) 



Figure 1: Gluino- mass dependence of the colour-decomposed NLO hard matching coefficients for 
the qq initiated channel (a) and the gg initiated channel (b). The squark mass has been set to 
TTtq — 1.2 TeV, while the common renormalization and factorization scale has been set equal to the 
squark mass. The top quark mass is taken to be m t = 172.9 GeV. 

contribution to the cross section vanishes because it yields a p-wave contribution, which 
vanishes at threshold. 

These matching coefficients have been checked numerically using PROSPINO [14] and 
agree within the numerical accuracy of PROSPINO. They also agree with the a\ coefficients 
presented in Ref. [35] to the percent level. 

4. Numerical results 

In this section we present numerical results for the NNLL-resummed cross sections with 
and without the Coulomb contributions. We show the results for squark-antisquark pair- 
production at the LHC for centre-of-mass energies of 7 TeV and 14 TeV. In order to evaluate 
hadronic cross sections we use the 2008 NLO MSTW parton distribution functions [55] with 
the corresponding a s (M§) = 0.120. We have used a top quark mass of m t = 172.9 GeV [56]. 
The numerical results have been obtained with two independent computer codes. 

It should be noted that the Coulomb effects can be screened by the width of the sparti- 
cles depending on the specific SUSY scenario. For consistency we will stick to the approach 
adopted in the NLO calculations, where this screening is not taken into account. In or- 
der to study the effects from the hard matching coefficients and the Coulomb corrections 
separately, we will compare several cross sections with the NLO result and discuss their 
contribution: 

• The NLL matched cross section is based on the calculations presented in [30-32] and 
will be denoted as <j NL °+ NLL . 

• The NNLL matched cross section without Coulomb contributions to the resumma- 
tion fj NL 0+ NNLL w/o Coulomb con t a i ns th e soft-gluon resummation to NNLL accuracy 
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matched to the full NLO result. The matching is performed according to Eq. (| 
The Coulomb correction to the resummation is not included, so C Coul >W in Eq. (| 
is set to zero. 

• The NNLL matched cross section (j nl O+nnll ^ oes i nc l uc l e the Coulomb contribution 
£>Coui,(i) f rom equation (|2.6| ). Also in this case Eq. (2^) has been used to match the 
cross section to the complete NLO result. 

The NLO cross sections are calculated using the publicly available PROSPIMO code [57], 
based on the calculations presented in Ref. [14]. The QCD coupling a s and the parton 
distribution functions at NLO are defined in the MS scheme with five active flavours. The 
masses of squarks and gluinos are renormalized in the on-shell scheme, and the SUSY 
particles are decoupled from the running of a s and the parton distribution functions. No 
top-squark final states are considered. We sum over squarks with both chiralities (qi 
and qa), which are taken as mass degenerate. The renormalization and factorization scales 
H are taken to be equal. 

We first discuss the scale dependence of the cross sections. Figure || shows the squark- 
antisquark cross section for rriq = nig = 1.2 TeV as a function of the renormalization and 
factorization scale \i. The value of \x is varied around the central scale [i§ = rriq from 
fi = /io/5 up to fj, = 5 no and results are shown for the LHC at CM energies of 7 TeV (a) 
and 14 TeV (b). 



0.0020 
0.0018 
0.0016 
0.0014 
0.0012 
0.0010 
0.0008 
0.0006 
0.0004 
0.0002 




a(pp^qq + X)[pb] 
VS = 7TeV 

fj,0 = rriq = nig = 1200 GeV 




LO •---..„ 
NLO 

NLO + NLL 

NLO + NNLL w/o Coulomb 
NLO + NNLL 



O.r, 



1 



a{pp^t qq + X) [pb] 
yfS = 14 TeV 



Ho = rriq = m g = 1200 GeV 



(a) 




Figure 2: The scale dependence of the LO, NLO, NLO+NLL and NLO+NNLL (both with the 
Coulomb part C Coul: W anc j without it) squark-antisquark cross sections for the LHC at 7 TeV (a) 
and 14 TeV (b). The squark and gluino masses have been set to rriq = m~ g = 1.2 TeV. 

For both collider energies we see the usual scale reduction going from LO to NLO. In- 
cluding the NLL correction and the NNLL contribution without the Coulomb part C Coul '^ 
improves the behaviour for moderate values of ///Mo> but a fairly strong scale dependence 
for small values of ///^o remains. Upon inclusion of the Coulomb corrections C 00 " 1 ^ 1 ) the 
scale dependence stabilises over the whole range. 
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Figure || shows the mass dependence of the scale uncertainty for the NLO, NLO+NLL 
and NLO+NNLL cross sections at the LHC. The squark and gluino mass have been taken 
equal and the scale has been varied in the range rriq/2 < fj, < 2rriq. As was to be expected 
from figure ||, the scale uncertainty reduces as the accuracy of the predictions increases. In 
the range of squark masses considered here, the NNLL resummation without the Coulomb 
corrections C Coul '^ already reduces the scale uncertainty to at most 10% for the LHC at a 
CM energy of 7 TeV and to even lower values for a CM energy of 14 TeV. The inclusion of 
the Coulomb term C Coul '^ in the resummed NNLL prediction results in a scale uncertainty 
of only a few percent for both collider energies. The effect of the threshold resummation is 
more pronounced for a collider energy of 7 TeV, which is due to the fact that the sparticles 
are produced closer to threshold in that case. 
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Figure 3: The scale uncertainty of the NLO, NLO+NLL and NLO+NNLL (both with the Coulomb 
part C Cou1 '^ and without it) squark-antisquark cross sections for the LHC at 7 TeV (a) and 
14 TeV (b). The common renormalization and factorization scale has been varied in the range 
w?/2 < fJ- < 2mq and the squark and gluino masses have been taken equal. 

Finally we study the K-factors with respect to the NLO cross section: 

K~ 



a 

NLO 



a 

where x can be NLO+NLL, NLO+NNLL w/o Coulomb or NLO+NNLL. In figure | we 
study the mass dependence of the K-factor for equal squark and gluino masses. 

At the central scale [i = rriq the JC-factor, and thus the theoretical prediction of the 
cross section, increases as more corrections are included. Also, the effect becomes more 
pronounced for higher masses. This was to be expected, since in that case the particles are 
produced closer to threshold. As can be seen in figure ||(a), the NNLL resummation without 
the Coulomb corrections C Cou ^^ already results in a 25% increase of the cross section with 
respect to the NLO cross section for squarks of 2 TeV and a CM energy of 7 TeV. The 
contribution from the Coulomb term to the resummed NLO+NNLL cross section is larger 
than the contributions provided by the 53 term in the exponential and the hard matching 
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Figure 4: The if-factor with respect to the NLO cross section of the NLO+NLL and NLO+NNLL 
(both with the Coulomb part C Cou1 '^ and without it) squark-antisquark cross sections for the LHC 
at 7 TeV (a) and 14 TeV (b) . The squark and gluino masses have been taken equal and the common 
renormalization and factorization scale has been set equal to the squark mass. 
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Figure 5: The i^-factor with respect to the NLO cross section of the NLO+NNLL squark- 
antisquark cross sections with and without the Coulomb contributions C 00 " 1 '' 1 ' for the LHC at 
7 TeV. Different ratios of the squark and gluino mass are shown, while the common renormaliza- 
tion and factorization scale has been set equal to the squark mass. 



coefficient yielding a total if-factor of 1.45. For the case of a CM energy of 14 TeV, 
which is shown in figure ||(b), the size of the NNLL contributions is smaller, since the 
sparticles are produced further away from threshold. However, for masses of 3 TeV, the 
iT-factor for the NNLL contribution without the Coulomb correction C Gonl '^ still yields 
1.13, whereas the inclusion of the Coulomb corrections increases this to 1.25. Although the 
effect from the Coulomb corrections could be somewhat smaller in reality due to the finite 
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lifetime of the squarks, figure [| suggests that the NNLL contribution will remain large. 

Figure || only contains the numbers for equal squark and gluino masses, but the effect 
of the gluino mass is small, as can be seen in figure |5[ In figure || the mass ratio r = mg/rriq 
has been varied. Although some effect can be seen, it is negligible compared to the size 
of the NNLL corrections. It turns out that this conclusion also holds for a collider energy 
of 14 TeV. Consequently the NNLL-resummed results are relatively independent of the 
relation between squark and gluino masses. 

The scale dependence of the cross section shows the best stability after including both 
the hard matching coefficients and the Coulomb contributions C Goui >^. This indicates 
that all these contributions should be taken into account to achieve the observed cancel- 
lation, see also figure 9 in [29]. However, the observed reduction in the scale dependence 
might be modified somewhat by the inclusion of the width of the particles or by match- 
ing to the full NNLO result, which is not available. In this context we note that, as a 
consequence of the NNLL accuracy of resummation, our matched cross section receives ad- 
ditional non-logarithmic NNLO contributions, which would have been consistently treated 
if matching to NNLO had been possible. A very conservative estimate of the scale un- 
certainty is provided by the NLO+NNLL w/o Coulomb results, which do not include the 
Coulomb corrections. 



5. Conclusions 

We have performed the NNLL resummation of threshold corrections for squark-antisquark 
hadroproduction. In particular, the previously unknown hard matching coefficient C^ 1 ' , 
needed at this level of accuracy, has been calculated analytically. We have also numerically 
evaluated the NNLL resummed cross section, matched to the NLO fixed-order expression, 
for squark-antisquark production at the LHC with CM energies of 7 and 14 TeV. At 
both collider energies the total cross section increases at the central scale. At 7 TeV 
collision energy and for a squark mass of 2 TeV, the NLO+NNLL squark-antisquark cross 
section is larger than the corresponding NLO cross section by as much as 45%. The 
correction is reduced to 25% if the contributions due to Coulombic interactions are not 
taken into account. In addition, the scale dependence is reduced significantly, particularly 
after inclusion of the Coulomb corrections. This information should be used to improve 
current limits on SUSY masses or, in the case that SUSY is found, to more accurately 
determine the masses of the sparticles. 
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A. The NNLL functions 



(3) 

In the following we list the explicit expression for the NNLL function g\/_^.= j, cf. Refs. [47, 
49]. Expressions for the LL and NLL functions can be found in Refs. [31,36]. 
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The NNLL functions q _ .= r and q 
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with A = boa s {ii 2 )\og{N), fj, the common renormalization and factorization scale, and je 
Euler's constant. The coefficients of the QCD beta function are denoted by b n and the first 
three coefficients are given by [59,60] 
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where rii denotes the number of light quark flavours, C4 = N c and Cp 



Nj-1 
2N, 



with N r the 



number of colours. The universal, process independent coefficients up to NNLL accuracy 
are given by [47] 



A ( 2) 



1 



(A.3) 



2 

(3) _ 1 



a 



Cli 2 ^- ? C (2) + ^C(3) + yC(2) 2 ) + C Fni (-g + 2C(3) 



24 9 ' w 6 

/ 209 10 7 A \ n, 

+ C *"> -108 + y« 2 »-3« 3 ) "27 



21 



and [42,61,62] 
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with the colour factor Cj = for i = q,q and Cj = Ca for i = g. The process dependent 
coefficients read [36] 
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B. Mellin transforms of the Coulomb corrections 

In this appendix we present the analytical results for the Mellin transforms of the Coulomb 
corrections in terms of the Euler beta function f3, the digamma function ^ and the hy- 
pergeometric functions 2-fi and 3F2. For the subprocess qiq~j — > qq the expressions for the 
colour-decomposed Coulomb part in iV-space are given by 
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whereas for the subprocess gg — > qq they read 
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Here r = m g /mq and n/ = 5 denotes the number of light quark flavours. The function 
ijv(r) is given by 
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4r 2 ' " (r 2 + l) 2 ' 

The Euler beta function f}(N + 1, 1/2) is abbreviated by 



(B.6) 



(B.7) 



B N = P(N + 1,1/2). 



C. How velocity factors can deform the phase space integration 

As mentioned in section |3] the dipole function given in Ref. [54] has been modified, which 
makes it unsuitable for our calculation. We will explicitly show the effect of this modifica- 
tion by considering the change in the second term of Eq. ( |3.2| ) , which is denoted as IgQ l Q 
in Eq. (5.23) of Ref. [54]. In Ref. [54] finite pieces of the integrand are taken into account 
as well, but since we showed in section || that the only contribution at threshold comes 
from the singular part of the integrand, we will omit these terms. 

The singular term of the integrand yields a 1/e-pole and a finite piece. The pole cancels 
the pole of the first term of the dipole function (|3.2| ), while the finite piece contributes to 
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the hard matching coefficient. In its unmodified form, the finite piece is given by 
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and 
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Pg-Pj +Pg-Pl+Pj - Pi l-2fi 2 

The approximation in the second line of Eq. (C.l) is suitable near threshold, where /i ~ 1/2. 
Exactly at threshold the finite part equals 4 — 41og(2) and exactly cancels the contribution 



from the first term of Eq. (3.2). 

In Ref. [54], the integrand in I^qq has been multiplied by velocity factors in order to 
simplify the integration: 



v gQ,Q 



(l-y)yi^V 



v gQ,Q V[V + (l-2^ 2 )(l-y)] 2 -4/x 2 y/(l-y/2)*-W 

where the approximation in the second step holds near threshold. The velocity factor 
effectively replaces \A — in the denominator of Eq. ( |C.1| ) by ^/(l — y/2) 2 — 4ju 2 , which 
amounts to a shift comparable in size to the value of the numerator. In the strict soft limit 
y vanishes and the velocity factors have no effect. However, we are integrating over gluons 
that are not soft compared to the energy above threshold y/s — 2m, so we also need the 
correct behaviour away from the strict soft limit. In fact, if the velocity factors are included 
the integral I^qq vanishes at threshold, so it does no longer cancel the contribution from 
the first term of Eq. (pDP . 

Usually the velocity factors do not pose a problem in calculations using dipole sub- 
traction, since the terms are subtracted from the real part and added to the virtual part. 
Therefore it does not matter if a dipole function is deformed, as long as the pole is re- 
produced. Finite contributions can always be moved between and a^ 3 \ However, 
we argued that ct^ 3 } vanishes at threshold due to phase-space suppression, which is not 
true if the phase space integration is deformed by velocity factors. Therefore we need the 
unchanged dipole function for this particular calculation. 

D. The hard matching coefficients for squark-antisquark production 

Here we present the exact expressions for the hard matching coefficients for the squark- 
antisquark production processes as defined in Eq. ([2.6;) . We sum over squarks with both 
chiralities (qi and (Jr). No top-squark final states are considered and all squarks are 
considered to be mass-degenerate with mass mx. Top squarks are taken into account in the 
loops, where they are taken to be mass-degenerate with the other squarks. The calculation 
is outlined in section || and was done with FORM [58]. We first define: 



4mim2 
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and 



2 2 2 
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2 2,2 

m + = m- Q + m~ 



where m g is the gluino mass. Denoting the number of light flavours by ni = 5, the total 
number of flavours by n/ = 6 and the number of colours by N c = 3, we also define: 
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We denote the factorization scale by hf, the renormalization scale by (ir and Euler's 
constant by 7^. For the qq — > qq process both the singlet 1 and the octet 8 representation 
contribute: 
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In this equation the last two lines are proportional to the quadratic Casimir invariant of 
the representation, which is zero for the singlet and iV c for the octet representation. We 
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have defined the functions: 
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For the — > qq process the antisymmetric octet 8a vanishes because it yields a rj-wave 
contribution, which vanishes at threshold. The hard matching coefficients for the singlet 1 
and the symmetric octet 8s do contribute: 
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where in the second equation the representation / can be the 1 or the 8s and the last line 
is proportional to the quadratic Casimir invariant of the representation. 
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